Geometric phases for wave packets in a uniform magnetic fieldQ 
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A wave packet of a charged particle always make cyclic circular motion in a 
uniform magnetic field, just like a classical particle. The nonadiabatic geomet- 
ric phase for an arbitrary wave packet can be expressed in terms of the mean 
value of a number operator. For a large class of wave packets, the geometric 
phase is proportional to the magnetic flux encircled by the orbit of the wave 
packet. For more general wave packets, however, the geometric phase contains 
an extra term. 
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I. INTRODUCTION 



Since the discovery of the geometric phase |T|, |], 0, the subject has been studied for various 
systems [Q, [7], |S| , especially particles with spin in time-dependent magnetic fields || [10], 
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TT[ |12| , [13| , |14| , |I5[ [T6| , [17| , p~8| , [0| . To our knowledge the geometric phases for wave packets 
that undergo cyclic motion in the configuration space were not considered in the literature. Of 
special interest is the wave packets of a charged particle moving in a uniform magnetic field, 
because this represents a rather simple and realizable physical situation. It should be noted 
that a geometric phase can occur even though the system is time independent, just like spin 
in a uniform magnetic field [0 . 

In this paper we consider a charged particle moving in a uniform magnetic field. In Sec. II 
we briefly review the energy levels and the eigenstates of the system. Though the solutions to 



the system is well known in the textbooks p0| , gj] , we present the results somewhat differently. 
In Sec. [HI] we consider the time evolution of an arbitrary wave packet. The time evolution is 
always cyclic and the center of the wave packet makes circular motion like a classical particle 
regardless of the initial condition. The nonadiabatic geometric phase in a cycle can be expressed 
in terms of the mean value of a number operator. This is similar to the case for a particle with 
spin in a rotating magnetic field, where the nonadiabatic geometric phase can be expressed in 
terms of the mean value of the component of the angular momentum along the rotating axis 



17 , lSfl . We also establish a linear relation between the nonadiabatic geometric phase and 



the magnetic flux encircled by the circular orbit of the wave packet. There is an inhomogeneous 
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term in this relation, which we call an extra term. A large class of wave packet is presented 
in Sec. [IV| for which the extra term vanishes and thus the nonadiabatic geometric phase is 
proportional to the magnetic flux. For wave packets that do not belong to this class, the extra 
term depends on the initial condition. Again this is similar to the case for particles with spin 
in a rotating magnetic field A brief summary is given in Sec. M. 



II. REVIEW OF THE MODEL 

Consider a charged particle with electric charge q and mass M, moving on the xy plane 
under the influence of a uniform magnetic field B = Be z where e z is the unit vector in the z 
direction and B is a constant which can be taken as positive without loss of generality. We 
take the symmetric gauge A x = —By/2, A y = Bx/2, then the Hamiltonian reads 

H = ^(P 2 X + Vl) + \muI{x> + y 2 ) - \eu B L z , (1) 

where uj b = \q\B/Mc > 0, L z = xp y — yp x , and e = 1 (—1) if q is positive (negative). We 
define 

Mu B x + i2p x Mu B y + i2p y 

di = , a 2 = ; 2 

2y/Mhu B 2y/MHu B 

The nonvanishing commutators among these operators and their Hermitian conjugates are 
fa,;, all = The Hamiltonian can be recast as 



H = -fkj B [(a\ai + a\a 2 + 1) + ie(a{a 2 — axal,)]. (3) 
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Next we define 



a = -^=(a! + iea 2 ), 6 = -^=(ai - iea 2 ), (4) 

and their Hermitian conjugates. They satisfy 

[o, a*] = [6,6+1=1, (5) 
and all other commutators vanish. In terms of these operators, we have 

H = hio B (a j a+ |), (6) 

and 

L z = eh{b ] b-a ] a). (7) 

We see that the Hamiltonian becomes that for a simple harmonic oscillator. Thus the 
energy levels are 



E n = tuo B {n+ |), n = 0,l,2, 



(8) 
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One can take the eigenstates to be the common ones of N = a^a and N' = tfb, denoted by 
\nn'), satisfying 

N\nn') — n\nn'), N'\nn') — n'\nn'), n, n' = 0, 1, 2, . . . . (9) 

They are also common eigenstates of H and L z , with eigenvalues E n and e(n'—n)h, respectively. 
These states are given by 

\nn') = -^=(at)»(6t)»'|00>, (10) 
v nln'l 

where the ground state 1 00) satisfies 

a|00) = 6|00> = 0. (11) 

One can work out the wave functions in the configuration space for these eigenstates and show 
that they are essentially the same as those obtained by solving the Schrodinger equation in 
the cylindrical coordinates. However, we are not interested in the quantum number n' in this 
paper, so we will only deal with the quantum number n. We consider the eigenstates \n) of N 
or H, satisfying 

N\n)=n\n), H\n)=E n \n), n = 0,1,2,.... (12) 
These are also called number states in the following. They are given by 

|n) = 4r( at ) n l°>' ( 13 ) 

Vnl 



where the ground state |0) satisfies 

a|0)=0. (14) 
Obviously, the state \n) is a linear combination of \nn'), that is 

oo 

J2 C n'\nn), (15) 



\n) = 

n'=0 
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where the coefficients C n > are arbitrary except satisfying ^2^ =0 \C n /\ = 1 such that \n) is 
normalized. Therefore there must exist some big freedom in the wave function for the state 
|n). 

We define a complex number z and its complex conjugate z* as 



z = ^ (x+iey) , z > = ^ (x - iev) , (16) 

then in the configuration space we have 

a=^=(z + d*.), <f = +=(z* - d,). (17) 
The wave function for the ground state |0) is obviously 

M^z*)=exp(-z*z)f(z), (18) 
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where the function f(z) is such that ipo(z, z*) is well behaved everywhere and is normalizable. 
Thus f(z) is rather arbitrary. In particular, any polynomial satisfies the requirement. We 
assume that ipo(z, z*) has been normalized, then the normalized wave function for the higher 
excited state \n) is 



1pn{z,Z*) 
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exp(zrz)dg[exp(-2z*z)f{z)]. 



(19) 



Therefore the freedom in these eigenstates lies in the arbitrariness of f(z). Of course this 
freedom corresponds to the degeneracy of the Landau energy levels (|8]). An arbitrary state of 
the system can be expressed as a linear combination of the above number states. 

A class of states (or wave packets in the configuration space) that are of special interest in 
the following are the so called displaced number states. Similar to those for a simple harmonic 
oscillator |2T], they are defined as 

\n, a) = D(a)\n), 
where D(a) is the unitary displacement operator 



(20) 



D(a) = exp(cW — a* a), 



(21) 



where a is a complex number. One can define a more general displaced state \(p, a) by acting 
D(a) on an arbitrary state \ip). An important property is that if the wave function for \ip) is 
ip(z,z*) or <f(x,y), then that for \<p,a) is 



fa(z,z*) = exp 



az 



a z 



V2 



if z 



a 
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(22a) 



or 



2h 



Muj 



-a x ,y — e 



B 



2h 



Muj B ay 



(22b) 



<Pa(x, y) = exp iy ~^-( a y x ~ ea xV) f \ ' - 

where a x = Re a, a y = Ima. Thus the wave packet of the displaced state is essentially an 
entire displacement of the original one, and the above definition for the displacement operator 
D(a) seems to be sound. 



III. THE GEOMETRIC PHASE FOR AN ARBITRARY WAVE PACKET 

Consider the time evolution of an arbitrary wave packet. The initial state at t — may be 
expressed as 

oo 

l^(0)> = $>». (23) 

n=0 

where the coefficients c n are arbitrary except satisfying J2^=o \ c n\ 2 — 1 such that \ip(0)) is 
normalized. Since the Hamiltonian is given by Eq. (|6]), the state at a latter time t is 

oo 

ty(t)) = e- iHt ^(0)) = J2 c n^M~^Bt(n + i)]|n). (24) 

n=0 



5 



At the time T = 2it/ujb, we have 

\ij(T)} = e---mO)}- (25) 
Therefore any state is cyclic, and the total phase change in a cycle is 

5 = -7T, mod 2tt, (26) 
which is independent of the initial condition. The expectation value of H in the state \tp(t)) is 

(H) t = m)\H\m) = (mmm} = + ±), ^ 

where the mean value (N) is evaluated in the initial state. The dynamic phase is 

(3 = -fT 1 I (H) t dt = -7T - 2n(N). (28) 
Jo 

This depends on the initial state. The nonadiabatic geometric phase is 

7 = 5 - (3 = 2n(N), mod 2vr. (29) 

It is proportional to the mean value of the number operator N. This is similar to the case for 
a particle with spin in a rotating magnetic field, where the nonadiabatic geometric phase can 
be expressed in terms of the mean value of the component of the angular momentum along the 
rotating axis || |17], [18]]. A number state has obviously a vanishing geometric phase (modulo 
2tt) since its time evolution is trivial. 

On the other hand, we consider the motion of the center of the wave packet. The position 
of it is characterized by the mean value of the coordinate variables x and y. We denote 

& = <V(0)|ahK0)), Pb = (^(0)161^(0)). (30) 

Using the relations 

e im/n ae -iHt/n = aexp( _i^) ? e im/H be -im/H = 6j (31) 
which can be easily verified, we have 

(V(t)|a|V(t)> = (3 a exp(-iu B t), (ip(t)\b\i/;(t)) = f3 b . (32) 

This leads to 



2h / 2h 

£ * ~ ^ ~Mu^P hx = V Mo7^'^ a ' COS ^ Bt ~ ar gA0> ( 33a ) 



2/z, / 2/i 

Vt + c Vm^^ = " £ fe IAI Sin( ^ _ arg/?a) ' (33b) 

where x t = (ip(t)\x\ip(t)), y t = {i ) if)\y\' l l ) it)) ^ and ^ = Re /3b, (3^ = Im/3 b . Therefore the 
center of the wave packet always makes a circular motion. The angular frequency is ujb- The 
motion is clockwise when q > and anticlockwise when q < 0. These are all the same as for 
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a classical particle. The radius of the circle depends only on (3 a , while fib only determines the 
position of the center of the circle. The magnetic flux go through the area encircled by the 
circular orbit is defined as positive (negative) if the motion is anticlockwise (clockwise). This 
turns out to be 

$ = -~2n\p a \ 2 = --2n(a^)(a) = --2vr|(a)| 2 , (34) 
q q q 



where the expectation values are evaluated in the initial state. Compared with Eq. ( P9|) we 
obtain 

7 = _££ + 27r(Aa) 2 , mod 2tt. (35) 

he 

where 

(Aa) 2 = (a+a) - (a t )(a) = (a f a> - |(a>| 2 = ((a - (a)) t (a - (a))). (36) 

Thus the nonadiabatic geometric phase contains two terms, the first is proportional to the 
magnetic flux encircled by the orbit of the center of the wave packet, the second is an "extra 
term" . In Sec. [TV] we will show that the extra term vanishes (modulo 27r) if the wave packet 
is initially a displaced number state. For a more general wave packet, however, the extra term 
depends on the initial condition. Again this is similar to the case for particles with spin in a 
rotating magnetic field |19|] . 



IV. THE GEOMETRIC PHASE FOR A DISPLACED NUMBER STATE 

Consider the special case where the initial state is a displaced number state 

|^(0)) = \n,a) = D(a)\n). (37) 

Using the relations 

D ] (a)aD(a) = a + a, D ] (a)a ] D(a) = a 1 " + a*, (38) 
which can be easily verified, we have 

(a) = (V»(0)|a|V(0)> = a, (at) = (^(0)1^(0)) = a*, (39) 

and 

(a f a) = <V^(0)|a t a|V(0)> = n + \a\ 2 . (40) 

Therefore in these states 

(Aa) 2 = n, (41) 

and 

a$ 

7 = --i— , mod 2tt. (42) 

he 
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This holds for all numbers n = 0, 1, 2, . . . , and for any complex number a. Once again this 
result is similar to the case for spin in a time-dependent magnetic field where for special 
initial conditions the nonadiabatic geometric phase is always proportional to the solid angle 
subtended by the trace of the spin vector. (For a rotating magnetic field this is well known. 
For an arbitrarily varying magnetic field, see [|i~3"| , [H| , |T9|j.) 

It should be noticed that the time evolution of a displaced number state is rather simple. 
Using relations similar to Eq. (|31~|), it is easy to show that 

e- iHt/h D(a)e iHt/h = D(a t ), (43) 

where 

a t = aexp(-iuBt). (44) 
Thus for the initial condition (|3T|), the time evolution is given by 

\4>(t)) = exp[-iu B t(n + \)]D{a t )\n) = exp[-iu B t(n + |)]|n, a t ). (45) 

We see that it remains to be a displaced number state at all latter times, except that the 
displacement parameter varies with time, and an overall phase factor is gained. Since the 
displacement only changes the position of the wave packet but not the shape of it, we conclude 
that a displaced number state keeps its shape unchanged while it is making circular motion. 

V. SUMMARY 

In this paper we have studied the time evolution and the associated nonadiabatic geometric 
phase for a wave packet of a charged particle moving in a uniform magnetic field. Any state of 
this system is cyclic, and the center of it always makes circular motion like a classical particle. 
The nonadiabatic geometric phase in a cycle can be expressed in terms of the mean value 
of a number operator. A linear relation between the nonadiabatic geometric phase and the 
magnetic flux encircled by the circular orbit of the wave packet is established. For wave packets 
that are initially displaced number states the nonadiabatic geometric phase is proportional to 
the magnetic flux. For more general wave packets it contains an extra term which depends 
on the initial condition. This shows that the nonadiabatic geometric phase is not always 
proportional to the geometric object (the magnetic flux in the present problem is a geometric 
object since it is proportional to the area through with it penetrates) involved in the cyclic 
motion, but may be a more general function of the latter. Similar results were encountered in 
the problem of particles with spin moving in a rotating magnetic field ]nj. We also showed 
that a displaced number state keeps its shape unchanged while it is making circular motion. 
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